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Motivation 

Finde die Iris im Bild 




Motivation: 

Finde die Iris im Bild 



max l 



(Daugman) Sucht folgendes: 



(r,*o,l/o)| G * 0)*^;/ 



27vr 



ds 



Auf Deutsch: Wir suchen das Maximum eines 
weichgezeichneten Gradienten uber bestimmte (normierte) 

Kurvenintegrale 



i — ► viele (neue?) Begriffe: 
Maximum, 'Weichzeichner', (partielle) Ableitung, Normierung, 

Integral 



Was ist ein Maximum? 




a) Intuitiver Zugang: Ein Buckel, 'Da ist es eben am groBten' 



b) Exakter Zugang: 

Eine auf D C R erklarte Funktion / hat in a e D ein (globales) 
Maximum, wenn f(x) < f(a) fur alle x e D 



Wieder ein neuer Begriff: Funktion 



Was ist eine Funktion? 



Programmierer: 

'Man wirft was rein und kriegt was raus' zusammen mit einem 

Codefragment der Art 

float quadrat (float x) { 
return x*x; 

} 



Mathematiker: 

eine Funktion auf einer Menge X ist eine Vorschrift f, die jedem 
Element x e X in eindeutiger Weise eine Zahl f (x) zuordnet. 



Funktionen (cont.) 



in eindeutiger Weise 



Auf gut Deutsch: wirft man ein und dasselbe x zweimal in die 
Funktion rein, kommt auch zweimal das selbe wieder raus 



Fragen zu Funktionen? 



(partielle) Ableitung 



Relevant fur Funktionen in mehreren Variablen. Wir haben erst 
mal nur eine - vor den partiellen kommen die handelsublichen ;) 



Dazu brauchen wir Differentialrechnung 



Differentialrechnung 



Eine Funktion f : I —> C auf einem Intervall I heiBt 
differenzierbar in xq g /, wenn der Grenzwert 

l im f( x )~f Oo) 

existiert. Dieser heiBt dann Ableitung Oder Differentia /quotient 
von / in xq. Die Funktion heiBt differenzierbar im Intervall I, 
wenn sie in jedem Punkt xq e I differenzierbar ist. 



Eine schone Definition. Aber was bedeutet das? Was ist dieses 

lim? 



Limiten 

lim steht fur Limes (Grenzwert) 
Folgengrenzwerte: Die Eulersche Zahl e (Beispiel) 



e = lim ( 1 + - 



n. 



n=l: (1 + 1)"= (l + i) =2 

n=5: (l + i) 5 ~ 2.49 

n = 1000: (l + iroo) ^2.72 



Limiten 




Steigung einer Geraden 
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k _ Ay _ v\ - yo 

Ax xi — xq 
andere Version der gleichen Sache: 

y = kx + d 



Sekantensteigung 



Eine Sekante ist eine Gerade, die mit einer Kurve zwei 
Schnittpunkte Pq und Pi gemeinsam hat 




Auch eine Sekante ist eine Gerade Wie vorher, nur mit 

Funktionswerten statt y-Werten 

, _ Ay _ /Oi) - /Qo) 

Z\x x\ — xq 



Tangentensteigung 



Fehlt nur noch der Limes ;) 

Kurve - 



Yy 




-Taogente 



S *< I'M ^V-^ 



xi — ► x (umbenennen - vollig egal wie man's benennt). Gehen 
mit x immer naher an xq (bilden den Limes) und bekommen 

l im f( x )~f Oo) 

X^XQ X _ X q 

(entspricht der Definition!) Ableitung von f(x) wird meist als 

f\x) bezeichnet 



Ein Beispiel 



Unsere Funktion: f(x) — oo j oo — 

+ Ax 



|jm f( x )~f Oo) _ |jm Oo + ^x) -(x ) 



X^XQ x _ X q X^XQ /\ x 

Xq + 2Axxq + Ax 2 — Xq 2Axxq + Ar 



+ 2Axxq + Ax 2 - x£ 
Mm — ^ = Mm 

X^XQ /\ X X^XQ /\ X 

= Mm 2xn + Ax = 2xn = 2x 



Regeln 



Funktion 


Ableitung 


Funktion 
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x n 


(n — 


C 





In x 


l 

X 


sin x 


COS X 




ae ax 


COS X 


-sin x 



Mehr Regeln 



U(x) + g(x)]' = f'(x) + g'(x) 
[c*f(x)]' = c*f'(x) 

Lf O) * g(x)Y = f(x)'g(x) + f(x)g(x)' (Produktregel) 
[44]' = f( x yn(x)-Kx) g (xy (Q U otientenregel) 

[/(<?<»)]' = /'(<?<») + g'{x) (Kettenregel) 



Beispiele: Linearitat 



/ 2 / 

[f( x ) H~ 9\ x )\ — [ x H~ sm(x)] = 2x + cos(x) 
[c * /(a)]' = [4x 3 ]' = 3 * 4 * x 2 = 12x 2 



Beispiele: Produktregel 



[f(x) * g(x)]' = f(x)'g(x) + f(x)g(x)' 
[/(x) * g(x)] = [x * sm(x)] = 3x * sin(x) + x * cos(x) 



Beispiele: Quotientenregel 



r f&l y — /(x) / gf(x)-/(x)gf(x) / 

f(x) y _ 3x , _ 3(x + 1) - 3x 
g(x) x + 1 (x + l) 2 



Beispiele: Kettenregel 



[/(sO))]' = f'(g(x))g'(x) 
[f(9(x))Y = [e 2x ]' = 2e 2 * 



partielle Ableitungen 



Bei Funktionen in mehreren Variablen hat man mehrere 

Moglichkeiten: z.B. /(r, <p) = rsin(ip) 

_d_ 

dx' 



Partielle Ableitung: §pf{r,(p) Oder ^f(r,<p) 



partielle Ableitungen (cont.) 



jfrf( r ,v) heiBt: leite nur nach r ab und lasse den Rest in Ruhe 

d d 

-^-f(r,<p) = —(rsin(<p)) = sin(tp) 
or or 



Zwischenstand 



max (r,x ,yo) 



Ga (r) 



* 



Weichzeichner? 



max (r,x 0l yo) 



dr J 27vr 

Gaussian Blur: eigentlich Faltung der aktuellen Funktion mit 

einer Gaussverteilung 
Theorie: Faltung zweier Funktionen via Produkt der 
fouriertransformierten Funktionen, danach Rucktrafo des 

Ergebnisses 



Sorry 



Heute keine Fouriertrafos ;) 

fit) = [°° F(u)e iujt duj 

J—oo 

F(oj) = f°° f(t)e iut dt 

J—oo 



Numerischer Zugang 



Wir mussen eigentlich nur glatten - Durchschnitt pro je 5 
Funktionswerten z.B. nehmen reicht aus 



for (i=start ; i<stop; i++) { 

newy=l/5 * (f (i-2)+f (i-l)+f (i)+f (i+l)+f (i+2) ) 

} 



Integralrechnung 



rnax (r,x ,y ) 



(j(7 \T) * — V : US 



drJ 27rr 

Hier: Ringintegral (Integral uber geschlossene Kurve) 
Man braucht nicht unbedingt ein Ringintegral. My humble 
opinion: Anschaulicher mit Flachenintegralen 
Davor: erst mal handelsubliche Integrale ;) 



Flache unter einer Kurve 




Suchen die Flache unter der Kurve zwischen x=0 und x=2 
Idee: wir zerschneiden die Flache in Rechtecke 



Flache unter einer Kurve 




N 

Yl ( a i) * ( Xn ~ x n-l) 
n=l 



Flache unter einer Kurve 



Fur TV —> oo (Limes!) 

rd N 

/ f(x)dx := Mm ( a i) * ( x n ~ x n-l) 

Jc N^oo h 

n=l 

Je kleiner die Rechtecke (also desto mehr), desto naher kommt 
unsere Summe dem wahren Wert - dem bestimmten Integral 



Was heiBt integrieren? 



Mit der Integration lost man das Umkehrproblem, aus der 
Ableitung f(x) die Funktion f(x) zu bestimmen - also zu einer 
Funktion f(x) die Stammfunktion f(x) zu finden 

/O) = / f(x)dx 



Regeln? 



Differenzieren funktioniert immer - integrieren nicht! Es gibt 

nicht immer eine Stammfunktion 



Funktion 


Stammfunktion 
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Exkurs: Koordinatensysteme 



Kartesische Koordinaten: x, y 
Polarkoordinaten: Dazu brauchen wir etwas Trigonometrie 



Winkelfunktionen 



sin 




Winkelfunktionen (cont.) 



Weg mit Grad, her mit Radiant: 360 deg = 27r rad 

irrad 



deg = 



180 



sin(0) = 0, sm(|) = 1, sin(ir) — 0, sin(J^-) — — 1 
cos(0) = 1, cos(|) = 0, cos(tt) = -1, cos(^) = 



Exkurs: Koordinatensysteme 
(cont.) 



Kartesische Koordinaten: x, y 
Polarkoordinaten: 

x = rcos((p) 
y = rsin((p) 

(Kugelkoordinaten, Zylinderkoordinaten, ...) 



Ein Tip 



Nicht versteifen auf Integrale als Flachen unter einer Kurve! 
Besser: Integration und Differentation als Einheit sehen 



Flachenintegrale 



Flacheninhalt in einer Variable hatten wir schon. x war das 

Funktionsargument 
Aber wie berechnet man den Flacheninhalt eines Kreises? Dazu 
brauchen wir mehr als eine Variable und somit ein 

Flachenintegral 



Exkurs: Baby-Fubini 



Satz von Fubini: Vertauschung der Reihenfolge der Integration 

J i^J f(x,y)dx ) Jdy = J i^J f(x,y)dy ] Jdx 



Blodes Beispiel: Rechteck 




nutzliches Beispiel: Kreis 




Mit kartesischen Oder (viel besser) Polarkoordinaten 



A = 



rr r27T rr 

/ / rdrd(f — I 27rrdr = 
JO JO Jo 



— r 2 7r 



Beispiel: Kreisring 




Exkurs: Normierung 




Endlich! 



Wir haben soweit alle Mathe die wir brauchen, um das Problem 

(endlich) anzupacken :) 
Nochmal zur Erinnerung: 

d rl(x,y) 



max 



O,xo,2/o) 



Ga (r) 



* 



ds 



dr J 27rr 

Das Ringintegral schenken wir uns - wir nehmen einfach einen 

'unendlich dunnen' Kreisring 



Suchs, Struppi 




Fur alle P= (x,y)\ Bilde den Durchschnitt der Helligkeitswerte 
(Graustufen) entlang Kreisringen mit Breite 1 Pixel von r = 

bis r — Tjxiax 



Suchs, Struppi 




Fur alle P= (x,y)\ Bilde den Durchschnitt der Helligkeitswerte 
(Graustufen) entlang Kreisringen mit Breite 1 Pixel von r = 

bis r — Tjxiax 



So siehts aus 



Auftragen der Durchschnittswerte der Helligkeit gegen r 




20 40 60 



Implement ierung 



C und MbSDL - fpc2.c 
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